Local Density of States in a Dirty Normal Metal Connected to a Superconductor 
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A superconductor in contact with a normal metal not only induces superconducting correlations, 
known as proximity effect, but also modifies the density of states at some distance from the interface. 
These modifications can be resolved experimentally in microstructured systems. We, therefore, study 
the local density of states N(E, x) of a superconductor - normal metal heterostructure. We find 
a suppression of N(E, x) at small energies, which persists to large distances. If the normal metal 
forms a thin layer of thickness L n , a minigap in the density of states appears which is of the order of 
the Thouless energy ~ hD/L 2 ^. A magnetic field suppresses the features. We find good agreement 
with recent experiments of Gueron et al. 



I. INTRODUCTION 

A normal metal in contact with a superconductor 
acquires partially superconducting properties. Super- 
conducting correlations, described by a finite value of 
the pair amplitude ( £C )V , t( £C ))' penetrate some dis- 
tance into the normal metal. This proximity effect 
has been studied since the advent of BCS theory (see 
Ref. [I] and references therein). Recently, progress in low- 
temperature and microfabrication tecknQlogy has rekin- 
dled the interest in these propertiesB~13. Interference 
effects in a ndirty normal metal increase the Andreev 
conductanceCm. The effect of the superconductor on 
the level statistics of a small normal grain has been 
investigated^. 

Whereas the order parameter penetrates into the nor- 
mal metal, the pair potential A(x) vanishes in the ideal 
metal without attractive interaction. Since A yields the 
gap in the single-particle spectrum of a bulk supercon- 
ductor, the question arises how the spectrum of the nor- 
mal metal is modified by the proximity to the supercon- 
ductor. Recently, this question has been investigated ex- 
perimentally by Gueron et alE. In their experiment, the 
local density of states of a dirty normal metal in contact 
with a superconductor was measured at different posi- 
tions and as a function of an applied magnetic field. 

In this paper, we evaluate the local density of states 
N(E, x) of a superconductor - normal metal heterostruc- 
ture with impurity scattering in a varietst-pf situations. 
We generalize earlier theoretical workEJlij by applying 
the quasiclassical Green's function formalism and by in- 
cluding the effect of a magnetic.-field. We compare with 
the experiment of Gueron et alB and find good qualita- 
tive agreement with the experimental data both in the 
cases with and without a magnetic field. 



II. THE MODEL 

In the following we will consider geometries as shown 
in Fig. |l|. The superconductor is characterized by a finite 
pairing interaction A and transition temperature T c > 0. 
In the normal metal we take A = T c = 0. Here we restrict 
ourselves to the dirty (diffusive) limit, £ 3> l c \, where 



£ = (D/2A) 1 / 2 is the superconducting coherence length 
at T = and l c \ is the elastic mean free path. The latter 
is related to the diffusion constant via D = ^vplei- 

The density of states (DOS) of this inhomogeneous sys- 
tem can be derived systematically within ±he quasiclassi- 
cal real-time Green's functions formalism^. In the dirty 
limit the equation of motion for the retarded Green's 
functions Ge, Fe readslia 
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The diagonal and off-diagonal parts of the matrix Green's 
function, Ge and Fe, obey the normalization condition 



Ft = 1 



(2) 



which suggests to parameterize them by a function 
9{E,x) via Fe — sin(0) and Ge = cos(0). Inelas- 
tic scattering processes are accounted for by the rate 
r in = l/2T in , while scattering processes from param- 
agnetic impurities are described by the spin-flip rate 
r s f = l/2r s f. At low temperatures the former is very 
small (rj n ~ 10 _3 A), and will be neglected in the follow- 
ing. 



(a) 



(b) 




/ 

A(y) > 


\z 


/f y 


t" 




















S 






X 


N 



FIG. 1. Geometries considered in this article, (a) A strictly 
one-dimensional geometry, (b) A more realistic geometry sim- 
ilar to experimental setup. 

For the geometry shown in Fig. [I] the order parameter 
can be taken real. On the other hand, in the vicinity 
of an N-S boundary the absolute value of the order pa- 
rameter is space dependent, and has to be determined 
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self-consistently. The self-consistency condition is con- 
veniently expressed in the imaginary-time formulation, 
where 

A(x)H-^-) = 2nT Fi»M ~ — ■ (3) 

Here, 10^ — 7rT(2/x -I- 1) are Matsubara frequencies. The 
summation is cut off at energies of the order of the Debye 
energy. The coupling constant in S has been eliminated 
in favor of T c , while the coupling constant in N is taken 
to be zero. 

In the case where the interface between N and SJaas 
no additional potential, the boundary conditions are£3 

F s (0_) =F E (Q+) (4) 

gSq7)^ Fe{q - ] = cffe) ^ Fe{0+) ■ 

Here, <J n ( s ) are the conductivities of the normal metal 
and the superconductor, respectively. The complete self- 
consistent problem requires a numerical solution. Start- 
ing from a step-like model for the order parameter, self- 
consistency was typically reached within 10 steps. Fi- 
nally the DOS is obtained from N(E) = N a ReG E (x), 
where No is the Fermi level DOS in the normal state. 
We will present now results for three different cases: 

A. The DOS near the boundary of a semi-infinite nor- 
mal metal and superconductor. 

B. The DOS in a thin normal film in contact with a 
bulk superconductor. 

C. The effect of a magnetic field on the DOS in an 
experimentally realized N-S heterostructure. 

In the following sections energies and scattering rates will 
be measured in units of the bulk energy gap A and dis- 
tances in units of the coherence length £ = (D/2A) 1 / 2 . 



Here 

oj = -iE + r in , 

d s = atan(— — — — ), 

-iE + r in 

. 0a -9. (-iE + T^) 1 / 2 . e 

sm -2~ = 7 (H£ + r in)2 + A 2 )V4 sm 2 ■ 

Several material parameters combine into the parameter 

7 = (oVi&M'&i) , (6) 

measuring the mismatch in the conductivities and the co- 
herence lengths of the two materials. Furthermore, £ s („) 
is defined by (D S („)/2A) 1 / 2 , where D s ( n ) is the diffusion 
constant of the superconductor (normal metal). 

The resulting DOS, N(E), in the normal metal at a 
distance x = 1.5£„ from the interface is shown in Fig. ^ 
for different values of the parameter 7. It shows a sub- 
gap structure with a peak below the superconducting gap 
energy E < A and a strong suppression at zero energy. 
The modification of the DOS is most pronounced at small 
values of 7 and at small distances. The smaller the en- 
ergy, the larger is the distance where the modifications 
are still visible. In particular at E — the DOS vanishes 
for all values of x. Pair-breaking effects lead to a finite 
zero-energy DOS, as will be shown later. 

1.5 | , 1 , 1 1 




III. RESULTS AND DISCUSSION 
A. DOS in an Infinite System 
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We assume that the normal metal and the super- 
conductor are much thicker than the coherence length 
L s ,L n ^> £ and investigate how the DOS changes con- 
tinuously from the BCS form 7V BC s(^)/^o = \E\ /{E 2 - 
A 2 ) 1 / 2 deep inside the superconductor to the constant 
value Nn(E)/No = 1 in the normal metal. 

In a first approximation, neglecting self-consistency 
and paramagnetic impurities, we can solve Eq. ^ ana- 
lytically, with the result 

4atan[ tan(0 o /4) exp(-y / 2ui/D n x) ] x > 
6(E,x) = { s +4atan[tan((0o-0 s )/4)x 



exp {J2Vuj 2 + A 2 /D s x) ] x<0. 

(5) 



FIG. 2. DOS in the normal metal at x = 1.5£„. 

Next we solve the problem self-consistently and present 
some numerical results for the case 7=1. We first con- 
centrate on the superconducting side of the boundary. 
As shown in Fig. | the peak in the DOS is strongly 
suppressed, changing from a singularity to a cusp, but 
it remains at the same position A as one approaches 
the boundary. On the other hand, the density of states 
with energies below A increases. The states with ener- 
gies well below A decay over a characteristic length scale 



£> s /(2VA 2 -£: 2 ), see Eq. (|). 
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FIG. 3. Density of states on the superconducting side of 
the N-S boundary. The inset shows the self-consistent pair 
amplitude. 

The DOS on the normal side at different distances from 
the NS-boundary is shown in Fig. |J. The pronounced 
sub-gap structure found in the approximate solution is 
still present in the self-consistent treatment. The fig- 
ure shows how the peak height and position change with 
the distance. In the absence of pair-breaking effects the 
DOS vanishes at the Fermi level for all distances (dotted 
curves). Inclusion of a pair- breaking mechanism (solid 
curves) regularizes the DOS at the Fermi level, and also 
the peak height is suppressed. The curves are in qualita- 
tive agreement with experimental data shown in Ref. ^. 
The self-consistent calculation presented here leads to a 
slightly better fit than the theoretical curves shown in 
Ref. |6| where a constant pair potential was used in the 
solution of the Usadel equation. In particular, the low- 
energy behavior of the experimental curves is reproduced 
correctly. 
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FIG. 4. Density of states on the normal side of an N-S 

boundary for two spin-flip scattering rates: F s { — (dotted 
lines) and r s f = 0.015 A (solid lines). 

At finite temperatures (but T <C T c ) we expect no 
qualitative changes in the behavior described above ex- 
cept that the structures in the DOS will be smeared out 
by inelastic scattering processes. Hence for an experi- 
mental verification temperatures as low as possible would 



be most favorable. 



B. DOS in thin N-layers 

Next we consider a thin normal layer in contact with 
a bulk superconductor, L s L n ~ £. The boundary 
condition at chosen to be d9(E,x)/dx = 0, 

i.e., the normal metal is bounded by an insulator. In this 
case the DOS on the N side develops a minigap at the 
Fermi energy, which is smaller than the superconducting 
gap. If the thickness of the normal layer is increased, the 
size of this minigap decreases. Results obtained from the 
self-consistent treatment are shown in Fig. ||. Details of 
the shape of the DOS depend on the location in the N- 
layeiO. However, the magnitude of the minigap is space- 
independent, as shown in the inset of Fig. ||. The magni- 
tude of the gap is expected to be related to the Thouless 
energy D/L^, which is the only relevant quantity which 
has the correct dimension. Of course the relation has to 
be modified in the limit L n — > 0. Indeed as shown in 
Fig. [| a relation of the form E g ~ (const £ + L n )~ 2 fits 
quite well. The sum of the lengths may be interpreted as 
an effective thickness of the N-layer since the quasiparti- 
cle states penetrate into the superconductor to distances 
of the order of £. The effect of spin-flip scattering in the 
normal metal on the minigap structure is also shown in 
the inset of Fig. [|. The minigap is suppressed as T s f is in- 
creased until a gapless situation is reached at T s { « 0.4 A. 
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FIG. 5. Minigap E g as a function of the normal-layer thick- 
ness. Inset: local DOS of an N-layer of thickness L n = 1.1£ 
in proximity with an bulk superconductor. 

We would like to mention that— a similar feature had 
been found before by McMillanEl within a tunneling 
model ignoring the spatial dependence of the pair am- 
plitude. We have considered here the opposite limit, 
assuming perfect transparency of the interface but ac- 
counting for the spatial dependence of the Green's func- 
tions. For r s f = our results for the structure of the 
DOS agree further with previous findings of Golubov and 
Kupriyanovuj and Golubov et al$3. Recently, a minigap 
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in a two-dimensional electron gas. in contact to a super- 
conductor has also been studiedlia. 



C. Density of states in a magnetic field 

An applied magnetic field suppresses the superconduc- 
tivity in both superconductor and normal metal. To 
study the effect of the magnetic field on our system we 
consider the geometry shown in Fig. |l|b. Because in the 
experimental setup the thickness of the films is much 
smaller than the London penetration depth, we can ne- 
glect the magnetic field produced by screening currents. 
Therefore it is reasonable to assume a constant magnetic 
field, which is present in both S and N. The vector po- 
tential is then chosen to be 



A = A(y)e x 



A(y) = Hy 



(7) 



Eq. ([!]) can be considerably simplified in the case that 
the size of the system in y-direction is smaller or of the 
order of £. The system is limited to — W/2 < y < W/2, 
where W ~ £. Therefore the Green's functions do not 
depend on y and the equation can be averaged over the 
width W. The equation reduces to the effective one- 
dimensional equation 



— (Ge 9 2 Fe — Fe 9^.Ge) — 

i-iE + T in ) F E - A G E + 2r cff G E F E 



(8) 



Here, r c ff = r„f + De 2 H 2 W 2 / 12 acts as an effective pair- 
breaking rate, which depends on the transverse dimen- 
sion and the applied magnetic field. 

If we approximate the Green's functions in the super- 
conductor by their bulk values, the DOS in the normal 
metal at zero energy can be calculated analytically: 



N(0) 
N 



tanh(2 v / r cff /£) x) 2r cff < A 
(l-a 2 )/(l + a 2 ) 2T cS >A 



where 



a 



Aexp(-2 v /r cff /£> x) 

2r eff + V 4r 2 ff — A 2 



(9) 



(10) 



In Fig. ^ the dependence of the DOS on r c ff at x = 
1.5£ is shown for two different spin- flip scattering rates 
(equal rates for normal metal and superconductor). At 
r c fr = 0.5A the field dependence of the DOS shows a 
kink. This kink arises because above this value of r e ff 
the zero-energy DOS in the superconductor is nonzero 
(gapless behavior), which leads to an even stronger sup- 
pression of the proximity effect. 
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FIG. 6. Zero-energy DOS in the normal metal at x = 1.5£ 
as a function of r e g. 

Figure [7] shows a quantitative comparison of these pa- 
suits with experimental data taken by the Saclay groupliil 
In this experiments, the differential conductance of three 
tunnel junctions attached to the normal metal part of 
the system was used to probe the DOS at different dis- 
tances from the superconductor. Accordingly, we have 
calculated the self-consistent DOS in the presence of a 
magnetic field throughout the system for paJl energies 
and determined the differential conductancecd. We used 
x = 1.8£, consistent with an estimate from a SEM- 
photograph, and used a spin-flip scattering rate of r s f = 
0.015A in the normal metal as a fit parameter. This is 
necessary in the framework of our approach to explain 
the finite zero-bias conductance at zero field. We, fur- 
thermore, assumed ideal boundary conditions at the NS 
interface, i.e., 7 = 1, the motivation being that great care 
was used in the experiment to produce a good metallic 
junction, and significant Fermi velocity mismatches are 
not to be expected. 

At low and high voltages the agreement with the ex- 
perimental data is good for all three field values. On the 
other hand, the maximum in the DOS is not reproduced 
well by our calculation. Including the effect of a non- 
ideal boundary, i.e., 7 < 1 leads to an increase of the 
peak in the DOS but to a less satisfactory fit at low volt- 
ages. We cannot resolve this discrepancy, but we would 
like to point out that our theory is comparatively sim- 
ple and does not include all the geometric details of the 
experiment (e.g., the geometry of the overlap junction 
is not really one-dimensional and would be difficult to 
treat realistically). Our intention is to show that the- 
oretical treatment described here contains the physical 
ingredients to explain the basic features of the experi- 
mental data. The overall agreement between theory and 
experiment demonstrated in Fig. 7 shows this to be the 
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FIG. 7. Quantitative comparison of experimenlH 

(dotted lines) and theory (solid lines). The experi- 
mental magnetic fields are H = 0, 400, and 800G; 
h = HeW(D/12A) 1/2 . The theoretical curves have been nor- 
malized by R, = dI/dV cxp {eV/A = 1.5). 



IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we have given a theoretical answer to the 
question asked in the introduction, viz., what is - beyond 
the proximity effect - the effect of a superconductor on 
the spectrum of a normal metal coupled to it. Using the 
(real-time) Usadel equations, we have calculated the lo- 
cal density of states in the vicinity of an N-S boundary in 
both finite and infinite geometries. It shows an interest- 
ing sub-gap structure: if the normal metal is infinite, the 
density of states is suppressed close to the Fermi energy, 
but there is no gap in the spectrum. This is the behavior 
found in a recent experiments. In thin normal metals we 
find a mini-gap in the density of states which is of the 
order of the Thouless energy. We have also investigated 
the suppression of these effects by an applied magnetic 
field and find good agreement with experiment. 
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